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The existence of the five-quark Fock states for the intrinsic charm quark in the nucleons was suggested some 
time ago, but conclusive evidence is still lacking. We generalize the previous theoretical approach to the light- 
quark sector and study possible experimental signatures for such five-quark states. In particular, we compare the 
d—u and u + d — s — s data with the calculations based on the five-quark Fock states. The qualitative agreement 
between the data and the calculations is inteipreted as evidence for the existence of the intrinsic light-quark sea 
in the nucleons. The probabilities for the \uuduu) and \uuddd) Fock states are also extracted. 
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The possible existence of a significant uudcc five-quark 
Fock component in the proton was proposed some time ago 
by Brodsky, Hoyer, Peterson, and Sakai (BHPS) yj to explain 
the unexpectedly large production rates of charmed hadrons at 
large forward x p region. In the light-cone Fock space frame- 
work, the probability distribution of the momentum fraction 
(Bjorken-x) for this nonperturbative "intrinsic" charm (IC) 
component was obtained [1]. The intrinsic charm originating 
from the five-quark Fock state is to be distinguished from the 
"extrinsic" charm produced in the splitting of gluons into cc 
pairs, which is well described by QCD. The extrinsic charm 
has a "sea-like" characteristics with large magnitude only at 
the small x region. In contrast, the intrinsic charm is "valence- 
like" with a distribution peaking at larger x. The presence of 
the intrinsic charm component can lead to a sizable charm pro- 
duction at the forward rapidity (xf) region. 

The x distribution of the intrinsic charm in the BHPS model 
was derived with some simplifying assumptions. Recently, 
Pumplin ifj showed that a variety of light-cone models in 
which these assumptions are removed would still predict the x 
distributions of the intrinsic charm similar to that of the BHPS 
model. The CTEQ collaboration y)] has also examined all rel- 
evant hard-scattering data sensitive to the presence of the IC 
component, and concluded that the existing data are consis- 
tent with a wide range of the IC magnitude, from null to 2-3 
times larger than the estimate by the BHPS model. This result 
shows that the experimental data are not yet sufficiently ac- 
curate to determine the magnitude or the x distribution of the 
IC. 

In an attempt to further study the role of five-quark Fock 
states for intrinsic quark distributions in the nucleons, we have 
extended the BHPS model to the light quark sector and com- 
pared the predictions with the experimental data. The BHPS 
model predicts the probability for the uudQQ five-quark Fock 
state to be approximately proportional to 1/ rriQ, where itiq is 
the mass of the quark Q [1]. Therefore, the light five-quark 
states uuduu and uuddd are expected to have significantly 
larger probabilities than the uudcc state. This suggests that 
the light quark sector could potentially provide more clear ev- 
idence for the roles of the five-quark Fock states, allowing the 
specific predictions of the BHPS model, such as the shape of 



the quark x distributions originating from the five-quark con- 
figuration, to be tested. 

To compare the experimental data with the prediction based 
on the intrinsic five-quark Fock state, it is essential to separate 
the contributions of the intrinsic quark and the extrinsic one. 
Fortunately, there exist some experimental observables which 
are free from the contributions of the extrinsic quarks. As 
discussed later, the d — u and the u + d — s — s are examples 
of quantities independent of the contributions from extrinsic 
quarks. The x distribution of d — u has been measured in a 
Drell-Yan experiment [3]. A recent measurement of s + s in a 
semi-inclusive deep-inelastic scattering (DIS) experiment 
also allowed the determination of the x distribution of u + 
d — s — s. In this paper, we compare these data with the 
calculations based on the intrinsic five-quark Fock states. The 
qualitative agreement between the data and the calculations 
provides evidence for the existence of the intrinsic light-quark 
sea in the nucleons. 

For a \uudQQ) proton Fock state, the probability for quark 
i to carry a momentum fraction Xi is given in the BHPS 
model yj] as 
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where the delta function ensures momentum conservation. N$ 
is the normalization factor for five-quark Fock state, and mi 



is the mass of quark i. In the limit of 1114^ >> 
where m p is the proton mass, Eq. [Tjbecomes 
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where N$ = N^/mf 5 . Eq. [2] can be readily integrated over 
x\, X2, X3 and £4, and the heavy-quark x distribution is: 



P{x 5 )= l -N b xl[\(l 



x 5 )(l + lOirs + x\) 
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One can integrate Eq. [3] over x$ and obtain the result V£ c = 
./V5/36OO, where is the probability for the \uudcc) five- 
quark Fock state. An estimate of the magnitude of V§ c was 
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FIG. 1: The x distributions of the intrinsic Q in the uudQQ config- 
uration of the proton from the BHPS model [1]. The solid curve is 
plotted using the expression in Eq. [5] for c. The other three curves, 
corresponding to c, s, and d in the five-quark configurations, are 
obtained by solving Eq. [T] numerically. The same probability V®® 
(jyQQ _ o.Ol) is used for the three different five-quark states. 



given by Brodsky et al. yj as « 0.01, based on diffractive 
production of A c . This value is consistent with a bag-model 
estimate [5[]. 

The solid curve in Fig. Q] shows the x distribution for the 
charm quark (P(x§)) using Eq. [3] assuming V§ c = 0.01. 
Since this analytical expression was obtained for the limiting 
case of infinite charm-quark mass, it is of interest to com- 
pare this result with calculations without such an assumption. 
To this end, we have developed the algorithm to calculate 
the quark distributions using Eq. Q] with Monte-Carlo tech- 
niques. The five-quark configuration of {x±, ...,£5} satisfy- 
ing the constraint of Eq. Q] is randomly sampled. The prob- 
ability distribution P{xi) can be obtained numerically with 
an accumulation of sufficient statistics. We first verified that 
the Monte-Carlo calculations in the limit of very heavy charm 
quarks reproduce the analytical result for P(x<=,) in Eq. [3] 
We then calculated P(x§) using m u — md — 0.3 GeV/c 2 , 
m c = 1.5 GeV/c 2 , and m p = 0.938 GeV/c 2 , and the re- 
sult is shown as the dashed curve in Fig. Q] The similarity 
between the solid and dashed curves shows that the assump- 
tion adopted for deriving Eq.|3]is adequate. It is important to 
note that the Monte-Carlo technique allows us to calculate the 
quark x distributions for other five-quark configurations when 
Q is the lighter u, d, or s quark, for which one could no longer 
assume a large mass. 

As mentioned above, the insufficient accuracy of existing 
data as well as the inherently small probability for intrinsic 
charm due to the large charm-quark mass make it difficult to 
confirm the existence of the intrinsic charm component in the 
proton. On the other hand the five-quark states involving only 
lighter quarks, such as \uuduu), \uuddd), and \uudss), might 
be more easily observed experimentally. We have calculated 



the x distributions of the s and d quarks in the BHPS model for 
the \uudss) and \uuddd) configurations, respectively, using 
Eq.Q] The mass of the strange quark is chosen as 0.5 GeV/c 2 . 
In Fig. Q] we show the x distributions of s and d, together 
with that of c. In order to focus on the different shapes of the 
x distributions, the same value of vf® is assumed for these 
different five-quark states. Figure [T]shows that the x distribu- 
tions of the intrinsic Q shift progressively to lower x region as 
the mass of the quark Q decreases. The x distributions of Q 
originating from the gluon splitting into quark-antiquark pair 
(.9 QQ) QCD processes are localized at the low-x region. 
FigureQ]illustrates an important advantage for identifying the 
IC component, namely, the intrinsic charm component is bet- 
ter separated from the extrinsic charm component as a result 
of their different x distributions. Nevertheless, the probabil- 
ity for intrinsic lighter quarks are expected to be significantly 
larger than for the heavier charm quark. The challenge is to 
identify proper experimental observables which allow a clear 
separation of the intrinsic light quark component from the ex- 
trinsic QCD component. As we discuss next, the quantities 
d(x) — u(x) and u(x) + d(x) — s(x) — s(x) are suitable for 
studying the intrinsic light-quark components of the proton. 

The first evidence for an asymmetric u and d distribu- 
tion came from the observation H that the Gottfried Sum 
Rule |01 was violated. The striking difference between the 
d and u distributions was clearly observed subsequently in the 
proton-induced Drell-Yan J3J |8|] and semi-inclusive DIS ex- 
periments 01 . This large flavor asymmetry was in qualitative 
agreement with the meson cloud model which incorporates 
chiral symmetry II 1 CHI . Reviews on this subject can be found 
inRefs. EHil- 

The d(x) — u(x) data from the Fermilab E866 Drell-Yan 
experiment at the Q 2 scale of 54 GeV 2 Q3D is shown in Fig. [2] 




FIG. 2: Comparison of the d(x) — u(x) data with the calculations 
based on the BHPS model. The dashed curve corresponds to the 
calculation using Eq. [7] and Eq. [5] and the solid and dotted curves 
are obtained by evolving the BHPS result to Q 2 = 54.0 GeV 2 using 
jj, = 0.5 GeV and fi — 0.3 GeV, respectively. 
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The d(x) —u(x) distribution is of particular interest for testing 
the intrinsic light-quark contents in the proton, since the per- 
turbative g — > QQ processes are expected to generate uu and 
dd pairs with equal probabilities and thus have no contribution 
to this quantity. In the BHPS model, the u and J are predicted 
to have the same x dependence if m u = m,d- It is important to 
note that the probabilities of the \uuddd) and \uuduu) config- 
urations, V™ 1 and V$ d , are not known from the BHPS model, 
and remain to be determined from the experiments. Non- 
perturbative effects such as Pauli -blocking [14] could lead to 
different probabilities for the \uuddd) and \uuduu) configu- 
rations. Nevertheless the shape of the d(x) — u{x) distribu- 
tion shall be identical to those of d(x) and u(x) in the BHPS 
model. Moreover, the normalization of d(x) — u(x) is already 
known from the Fermilab E866 Drell-Yan experiment as 

[ (d(x) - u{x))dx = 0.118 ± 0.012. (4) 
Jo 

This allows us to compare the d(x) — u(x) data with the cal- 
culations from the BHPS model, since the above integral is 
simply equal to V dd — V™, i.e. 

/ (d{x) -u(x))dx = Vf-Vf 1 = 0.118 ±0.012. (5) 
Jo 

Figure |2] shows the calculation of the d(x) — u(x) distri- 
bution (dashed curve) from the BHPS model, together with 
the data. The x-dependence of the d(x) — u(x) data is not in 
good agreement with the calculation. It is important to note 
that the d(x) — u(x) data in Fig. [2] were obtained at a rather 
large Q 2 of 54 GeV 2 Hi]. In contrast, the relevant scale, /i 2 , 
for the five-quark Fock states is expected to be much lower, 
around the confinement scale. This suggests that the apparent 
discrepancy between the data and the BHPS model calcula- 
tion in Fig.|2]could be partially due to the scale dependence of 
d(x) — u(x). We adopt the value of n = 0.5 GeV, which was 
chosen by Gliick, Reya, and Vogt [15] in their attempt to gen- 
erate gluon and quark distributions in the so-called "dynami- 
cal approach" starting with only valence-like distributions at 
the initial /i 2 scale and relying on evolution to generate the 
distributions at higher Q 2 . We have evolved the predicted 
d(x) — u(x) distribution from Q\ = [i 2 — 0.25 GeV 2 to 
Q 2 = 54 GeV 2 . Since d(x) — u(x) is a flavor non-singlet 
parton distribution, its evolution from Qq to Q only depends 
on the values of d(x) — u(x) at Qq, and is independent of 
any other parton distributions. The solid curve in Fig. |2] cor- 
responds to d(x) — u(x) from the BHPS model evolved to 
Q 2 = 54 GeV 2 . Significantly improved agreement with the 
data is now obtained. This shows that the x-dependence of 
d(x) — u(x) is quite well described by the five-quark Fock 
states in the BHPS model provided that the Q 2 -evolution is 
taken into consideration. It is interesting to note that an ex- 
cellent fit to the data can be obtained if fi = 0.3 GeV is cho- 
sen (dotted curve in Fig.|2|i rather than the more conventional 
value of /i = 0.5 GeV. We have also found good agreement 
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FIG. 3: Comparison of the x(d(x) + u(x) — s(x) — s(x)) data with 
the calculations based on the BHPS model. The dashed curve corre- 
sponds to the calculation using Eq.Q] and the solid and dotted curves 
are obtained by evolving the BHPS result to Q 2 — 2.5 GeV 2 using 
(i = 0.5 GeV and /i = 0.3 GeV, respectively. 

between the HERMES d(x)-u(x) data at Q 2 = 2.3GeV 2 H 
with calculation using the BHPS model. 

We now consider the quantity u(x) + d(x) — s(x) — s(x). 
New measurements of charged kaon production in semi- 
inclusive DIS by the HERMES collaboration |4] allow the ex- 
traction of x(s(x) + s(x)) at Q 2 = 2.5 GeV 2 . Combining this 
result with the x(d(x) + u(x)) distributions determined by the 
CTEQ group (CTEQ6.6) UM, the quantity x(u(x) + d(x) - 
s(x) — s(x)) can be obtained and is shown in Fig. [3] This 
approach for determining x(u(x) + d(x) — s(x) — s(x)) is 
identical to that used by Chen, Cao, and Signal in their recent 
study [ 17] of strange quark sea in the meson-cloud model lfl8ll . 

An interesting property of u + d — s — s is that the con- 
tribution from the extrinsic sea vanishes, just like the case for 
d — u. Therefore, this quantity is only sensitive to the intrinsic 
sea and can be compared with the calculation of the intrinsic 
sea in the BHPS model. We have 

u(x) + d(x) — s(x) — s(x) = 
P uu {xu) + P d3 (x d ) - 2P sS (x s ), (6) 

where P ( 3 ( 3(xq) is the x-distribution for Q in the \uudQQ) 
Fock state. Although the shapes of the intrinsic u, d, s, s dis- 
tributions can be readily calculated from the BHPS model, the 
relative magnitude of the intrinsic strange sea versus intrinsic 
non-strange sea is unknown. We have adopted the assump- 
tion that the probability of the intrinsic sea is proportional to 
l/m 2 Q , as stated earlier. This implies that /{\{V% U + 
■pf)) = mi/m 2 « 0.36 for m a = 0.3 GeV/c 2 and 
m s = 0.5 GeV/c 2 . With this assumption, we can now com- 
pare the x(u(x) +d(x) — s(x)~ s(x)) data with the calculation 
using the BHPS model, shown as the dashed curve in Fig. [5] 
The prediction of the BHPS model is found to be shifted to 
larger x relative to the data. This apparent discrepancy could 
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again partially reflect the different scales of the theory and the 
data. Since u + d — s — s is a flavor non-singlet quantity, we 
can readily evolve the BHPS prediction to Q 2 = 2.5 GeV 2 
using Qo = H = 0.5 GeV and the result is shown as the 
solid curve in Fig. [3] Better agreement between the data and 
the calculation is achieved after the scale dependence is taken 
into account. It is interesting to note that a better fit to the data 
can again be obtained with \i = 0.3 GeV, shown as the dotted 
curve in Fig. [3] 

From the comparison between the data and the BHPS cal- 
culation using fi = 0.5 GeV in Fig. [3] one can determine the 
sum of the probabilities for the \uuduu) and \uuddd) configu- 
rations, Y,Tf (= Vf+P£ a ). We found that SP| S = 0.471. 
Together with Eq. [5] we have 



0.294. 



(7) 



It is remarkable that the d(x) — u(x) and the d{x) + u(x) — 
s(x) — s(x) data not only allow us to check the predicted 
x-dependence of the five-quark \uuduu) and \uuddd) Fock 
states, but also provide a determination of the probabilities 
for these two states. As expected, the extracted values for the 
five-quark Fock states probabilities in Eq.|7]depends on the as- 
sumption for the probability of the \uudss). For the limiting 
case of Vf = 0, we obtain Vf 1 = 0.097 and = 0.215, 
which reflect the range of uncertainty of the extracted values. 
It is interesting to note that values obtained in Eq. [7] are con- 
sistent with the 1 / m,Q assumption for the probability of the 
\uudQQ) Fock state. If one uses the bag model estimate of 
V$ c ~ 0.01 150, the 1/rriQ dependence would then imply that 
7^ to be ~ O.Ol^ml/m 2 ,) ~ 0.25, consistent with the results 
ofEq.0 

In conclusion, we have generalized the existing BHPS 
model to the light-quark sector and compared the calculation 
with the d—u and u+d— s — s data. The qualitative agreement 
between the data and the calculation provides strong supports 
for the existence of the intrinsic u and d quark sea and the ade- 
quacy of the BHPS model. This analysis also led to the deter- 
mination of the probabilities for the five-quark Fock states for 
the proton involving light quarks only. This result could guide 
future experimental searches for the intrinsic s and c quark 
sea. This analysis could also be readily extended to the hy- 



peron and meson sectors. The connection between the BHPS 
model and other multi-quark models 1 19, 20| should also be 
investigated. 
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